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We present a new numerical scheme for solving the initial value problem for 
quasiequilibrium binary neutron stars allowing for arbitrary spins. We construct 
sequences of circular-orbit binaries of varying separation, keeping the rest mass 
and circulation constant along each sequence. The spin angular frequency of the 
stars is shown to vary along the sequence, a result that can be derived analyt- 
ically in the PPN limit. This spin effect, in addition to leaving an imprint on 
the gravitational waveform emitted during binary inspiral, is measurable in the 
electromagnetic signal if one of the stars is a pulsar visible from Earth. 



We develop a new formalism for the hydrodynamical part of the initial 
value problem (I VP) that has, as its most important feature, the capability 
of providing solutions for binaries with stars with arbitrary spins. The 
gravitational part of the IVP (i.e.; the solution of the Hamiltonian and 
momentum constraints for a quasiequilibrium circular orbit) is solved using 
the Wilson-Mathews conformal "thin-sandwich" approach . This method 
is typically based on restricting the spatial metric tensor to a conformally 
flat form and imposing a helical Killing vector to the spacetime to enforce 
the quasiequilibrium "circular orbit" condition. The details of the new 
formalism, as well as the numerical code tests, have been presented in 
Marronctti and Shapiro 2 . We construct sequences of quasiequilibrium 
orbits with constant rest mass and constant relativistic circulation along the 
stellar equator. The stars are assumed to be identical and obey an n = 1 
polytropic equation of state. These sequences, which mimic evolutionary 
inspirals outside the innermost stable circular orbit (ISCO) driven by (slow) 
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gravitational radiation emission, provide insight into the coupling of spin 
and orbital angular momenta in general relativistic binaries. 
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Figure 1. Ratio of spin vs orbital angular velocities as a function of the coordinate 
separation. The filled circles (open squares) correspond to the zero-circulation sequence 
with stellar compaction (m/R)^ = 0.14 (0.19). The solid lines show the analytic PPN 
prediction (see Appendix F of Marronetti and Shapiro 2 for a derivation). 

Figure 1 shows the evolution of the stellar spin during the inspiral for 
two zero-circulation (irrotational) sequences. We note an increase of the 
spin that can reach up to 10% of the orbital angular velocity, before the 
binary reaches the ISCO. Figure 2 shows the evolution of the orbital angular 
velocity. The crosses mark the position of the ISCO as determined by the 
full general relativistic hydrodynamical simulations described in Marronetti 
et al. 3 and shown in Figure 3. We see that the ISCO is located prior to 
the termination of the quasiequilibrium sequence in both cases. Note also 
that our results agree very well with post-Newtonian predictions at large 
separations. This variation of the spin of the star throughout the inspiral 
is a potentially observable electromagnetic effect when one of the neutron 
stars is a pulsar. This result takes on added significance in light of the recent 
discovery 5 of the highly relativistic double-pulsar system J0737-3039, in 
which one of the pulsars is spinning rapidly with a period of 23 ms. 

References 

1. J. R. Wilson and G. J. Mathews, Phys. Rev. Lett. 75, 4161 (1995). J. R. 
Wilson, G. J. Mathews, and P. Marronetti, Phys. Rev. D 54, 1317 (1996). 



3 



0.04 




Coordinate Separation 



Figure 2. Orbital angular velocities as a function of the coordinate separation for zero- 
circulation sequences. The caption is identical to Fig. 1. The crosses mark the ISCO 
as determined by the full general relativistic hydrodynamical simulations described in 
Marronctti et al. 3 . The solid lines correspond to a 2.5 PPN expansion, as detailed in 
Pati and Will 4 . 




Figure 3. Coordinate separation vs time for binaries with stellar compaction (m/R)^, = 
0.14 and zero-circulation that start at different separations. The curves labeled "Stable" 
(outside the ISCO) and "Merger" (inside the ISCO) bracket our estimation of the location 
of the ISCO. (From Marronctti et al. 3 .). Solid (dashed) curves represent high (low) 
resolutions and the filled circles mark the time of surface contact. 
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